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Abstract 

The scattering of a weakly bound (halo) projectile nucleus by a 
heavy target nucleus is investigated. A new approach, called the Un- 
correlated Scattering Approximation, is proposed. The main approx- 
imation involved is to neglect the correlation between the fragments 
of the projectile in the region where the interaction with the target 
is important. The formalism makes use of hyper-spherical harmonics, 
Raynal-Revay coefficients and momentum-localized wave functions to 
expand projectile channel wave functions in terms of products of the 
channel wave function of the individual fragments. Within this ap- 
proach, the kinetic energy and angular momentum of each fragment 
is conserved during the scattering process. The elastic, inelastic and 
break-up S-matrices are obtained as an analytic combination involv- 
ing the bound wave function of the projectile and the product of the 
S-matrices of the fragments. The approach is applied to describe the 
scattering of deuteron on ^^Ni at several energies. The results are 
compared with experimental data and continuum-discretized coupled- 
channels calculations. 
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1 Introduction 



In the last years one of the main interests in nuclear physics has been focused 
on the study of halo nuclei, i.e., weakly bound and spatially extended sys- 
tems where one or two particles (generally neutrons) have a high probability 
of being at distances larger than the typical nuclear radii (see refs. |[l|,|2| for 
a general review on these nuclei). The ability to produce secondary beams 
of halo nuclei opened new possibilities of investigating their structure. Two 
basic experimental probes involving high energy reactions have been devel- 
oped to study halo structure. The first one is to measure the momentum 
distributions of the fragments coming out after a collision with light stable 
nuclei The second probe treats the analysis of Coulomb break-up cross 

section when the nuclei are incident on highly charged targets 

The first type of reactions has been treated in detail in a series of pub- 
lications by the group of Aarhus jSHTll- Here, the simplest approach to 



understand halo nuclei fragmentation reactions involves the instantaneous 
removal of one of the particles from the few-body halo system. Within this 
approach, known as "sudden approximation", one assumes that the binding 
system is removed without disturbing the motion of the constituent particles. 
This approximation is only justified for reaction times much shorter than the 
characteristic time for the motion of the particles within the few-body sys- 
tem. The sudden approximation has been extensively applied to the study 
of three-body halo nuclei, and in particular to the Borromean systems, i.e., 
three-body systems where all two-particle subsystems are unbound 
Final interaction between the two non-disturbed spectators seems to play a 
crucial role in order to explain the narrow neutron momentum distributions 
measured. The participant-target interaction was first described considering 
only absorption. Further improvements have been included recently, treating 
the interaction between the target and each of the halo particles by means 
of a phenomenological optical potential ||15|. The total cross section is then 
obtained by adding the contributions from all the participants in the halo 
nucleus. Processes where two or three halo particles interact simultaneously 
with the target are neglected. This is consistent with the fact that the model 
is only accurate for the outer part of the wave function |11, 1^. This means 



that those geometric configurations where more than one halo particle get 
close to the target during the collision should be excluded. This shadowing 
effect has been treated in previous works under different approaches. In the 
analysis of the Aarhus group the shadowing is accounted for by excluding the 
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participant wave function inside spheres around the two spectators [|TT|,|T^ ■ 
The second type of probe to study the structure of halo nuclei is by means 
of Coulomb elastic break-up reactions with a projectile, composed by a core 
and valence neutrons, incident on highly charged targets. The Surrey group 
has studied in detail elastic scattering of halo nuclei from target within the 
"adiabatic" approach, i.e., the intrinsic motion is very slow compared to the 
scattering motion [p!^ -[T9|. Moreover, the interaction between the projectile 
and the target is described considering only the interaction between the core 
and the target. This requirement is relevant to Coulomb dominated processes 
when the core is charged and the valence particle is neutral. In the case that 
strong interactions dominate, the above requirement is most likely to be valid 
when the number of core nucleons greatly exceeds the number of valence 



nucleons ||T^. Within these approximations, the elastic differential cross 
section factorizes into two terms, the cross section for a point-like projectile 
scattered by the target, and a form factor that contains the effects of the 
projectile structure. The range of validity of the adiabatic approximation 
is also discussed in |18| concluding that for a pure strong interaction 
the adiabatic approach is justified for a given projectile-target system at 
sufficiently high energy. On the contrary, in the case in which the Coulomb 
interaction dominates the validity of the "adiabatic" approach is questionable 
at forward scattering angles. 

Apart from the approaches mentioned, several other models have been 
proposed in the literature, starting from the pioneering work of Bang and 



Pearson ||20|, including eikonal |[2T|^3|, semi-classical [|J,|25|, and mixed ap- 
proaches to describe direct and sequential break-up ||2| 



The paper is organized as follows. In section 2 we develop the formalism of 
the uncorrelated scattering approximation (USA). Here the basic assumption 
is to neglect the correlation between the fragments in the region where the 
interaction with the target is strong. In this situation the orbital angular 
momenta and kinetic energies of the fragments are conserved during the 
collision process, and this leads to an analytic expression for the S-matrix 
of the composite system in terms of the S-matrices of the fragments. In 
section 3 we present a preliminary application of the developed approach to 
the case of elastic and break-up deuteron scattering on ^^Ni. In section 4 the 
conclusions are presented. 
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2 The Uncorrelated Scattering Approximation 



In this section we introduce a new approach to describe the scattering of a 
weakly bound nucleus by a heavy target. As it will be shown later, the S- 
matrices to the bound and break-up states of the composite system are given 
in terms of the ground state wave function and the S-matrices corresponding 
to the interaction of the fragments with tlie target. 

The interaction of a composite particle with the target can be expressed 
as the sum of two terms. On one side, an average force acting on the centre 
of mass of the projectile, which makes the projectile to scatter but does not 
excite or break it. On the other side, tidal forces that make the projectile 
rotate, excite or break up. Then, when a composite particle scatters from a 
target there are two opposite effects: i) the interaction between the fragments 
tending to keep the fragments bound, and ii) the tidal forces tending to break 
the system. In order to simplify our discussion, we assume that the mass of 
the target is much larger than the masses of the fragments. The Hamiltonian 
can be written then as 



p2 ^ 

H = — + ^ + VAB{r) + vat{Rat) + vbt{Rbt) (1) 

= ^ + ^ + VAB{r)+VAT{RAT)+VBT{RBT) (2) 

ZrriA 2,mB 

where M = tua + tub and m = mATriB/ {niA + mB). In this model, the 
interaction between projectile and target can be written as the sum of a fold- 
ing potential, vf{R) — {(pol^AT + vbt\<Po) , which does not affect the internal 
structure of the projectile, and a tidal potential, VT{R,r) = vat{Rat) + 
vbt{Rbt) — V£,'{R), which tends to break the projectile. The function |0o) 
describes the intrinsic ground state of the projectile. Note that for large 
distances R, the tidal forces, coming from the gradient of vt, are negligible 
compared to the force between the fragments, coming from the gradient of 
vab, that tend to keep them bound. Hence, it is reasonable to ignore the 
tidal forces for large distances. 

On the contrary, for small distances R, tidal forces can be large. In this 
reasonable approach is to ignore the force between the fragments, 
substituting the potential vab by a suitable constant v. For very tightly 
bound systems, tidal forces may not be strong enough to overcome the force 
between the fragments for any distance R. For these systems, scattering will 
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be predominantly elastic and governed by the folding potential. However, 
for weakly bound systems there will be a critical distance Rm below which 
tidal forces overcome the force between the fragments. The distance R„i can 
be associated to an angular momentum L^, so that Rm is the turning point 
of the wave- function, fulfilling 

^^g^ + = (3) 

Note that for L > Lm, tidal forces are not very important because the turn- 
ing point is beyond On the contrary, for L < L„i tidal forces will be 
important, and the correlation between the fragments may be neglected. 

Let us consider the situation in which tidal forces can be neglected. Thus, 
the Hamiltonian H, approximated by , can be decomposed as follows, 

= K + hR (4) 

hr = 7^ \-VAB{r) (5) 

hR = ^ + MR) ■ (6) 

The eigenstates of for a total energy E can be expanded in terms of prod- 
ucts of eigenstates of the internal Hamiltonian hr corresponding to energies 
e„, times eigenstates of Hr corresponding to energies E — en. We make use 
of a discrete and finite basis of N normalizable states of the relative motion 
of the fragments. These basis states include the bound states of the pro- 
jectile and the resonant states of the continuum. Diagonalizing the internal 
Hamiltonian hr in this basis, one obtains the eigenstates \nIM) with inter- 
nal energies e„. Thus, the energy of the relative motion of the projectile and 
target in the asymptotic region is En = E — Cn- The states that correspond 
to energies £"„ < do not contribute to the wave function asymptotically 
and, therefore, we restrict our basis space to En > 0. The states \nIM) are 
characterized by a given angular momentum /, M. Thus, we can write 

(r>/M) = <Pn{r)yiM{r) (7) 

(p>7M) = 4>nip)YiMip) (8) 
in coordinate and momentum space, respectively. 
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For the purpose of defining scattering magnitudes, let us consider the 
Hamiltonians free from the interaction with the target H^, where = + 
vf{R). The regular solutions of are characterized by a total energy E, 
orbital angular momentum L, internal angular momentum of the fragments 
/, total angular momentum J and internal energy e„: 

KiLjMjiE)) = MPnR)\n{LI)JMj), (9) 

where J^l represents a regular wave function of the free Hamiltonian that is 
just proportional to a Bessel function ji,(P„i?), whereas \n{LI)JMj) is the 
channel wave function in which the internal state with angular momentum / 
is coupled to the relative angular momentum L to produce the total angular 
momentum J, Mj. The momentum associated to the relative motion is given 
by P^/2M — E—En. If one constructs wave packets out of this wave function, 
one will have incoming waves for t ^ — cxo, and outgoing waves for t +00. 
The scattering in is such that the incoming waves will be unmodified, 
while the outgoing waves will be affected by the S-matrix due to the folding 
potential Sf{L, En), which will be diagonal in the channel basis. 

Let us now neglect the correlation between the projectile fragments. 
Then, the interaction vab is replaced by a constant v. The total Hamil- 
tonian can be written in terms of two non-interacting Hamiltonians 

H = hA + he + v (10) 

hA = 7r^ + VAT{RAT) (11) 

ZniA 

hB = :r^ + VBT{RBT) . (12) 

ZniB 

The eigenstates of H corresponding to an energy E can be expanded in terms 
of the product of eigenstates of /i^ and Kb, such that E — Ea+Eb+v. Given 
the adequate boundary conditions, it is straightforward to solve the scattering 
problem for the Hamiltonian H. We consider the Hamiltonian free from 
interactions with the target H'^, so that H = + vat{Rat) +vbt{Rbt)- A 
solution of this Hamiltonian is given by the product of regular wave functions 
in the co-ordinates Rat and Rbt, characterized by angular momenta La, Ma 
and Lb, Mb, and energies Ea and Eb. 

I*Lm.l.m,(^)) = JlAPaRat)JlAPbRbt)\LaMaLbMb). (13) 
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Then, constructing a wave-packet, we find that for t — > — oo, the wave func- 
tion is given by the product of incoming wave functions, while for t — +00 it 
is given by the product of outgoing wave functions. Cross terms containing 
the product of an incoming wave on one co-ordinate and an outgoing wave 
on the other are cancelled for t — > ±00. Then, if we switch on the inter- 
actions Vat + vbt, the incoming part is unaffected, while the outgoing part 
gets multiplied by the product of the elastic S-matrices generated by each 
potential. This means that the three-body S-matrix for the Hamiltonian H 
is diagonal in the basis characterized by the linear momenta and angular 
momenta of each fragment, and is given by the product of the S-matrices of 
each fragment S{La, Ea)S{Lb, Eb). 

A basic point in order to deal with the matching is to realize that the 
wave functions of can be characterized by the hyper-angular momentum 
K. In a basis of hyper-spherical harmonics, the wave functions obtained in 
the absence of interactions can be written as 

l^^KiLjAijiE)) = jK{Vn)\K{LI)JMj), (14) 

where Jk{x), that is proportional to the Bessel function Jk+2{x)^ is a regular 
solution of free three-body problem in terms of the hyper-radius 7^, given by 
TZ^ = + r'^m/M, and the hyper-momentum given by V'^/2M = E — v. 
The wave function \K{LI)JMj) can be written in terms of the hyper-angle 
a that defines the ratio of p to V, i.e., sin a = {^jM/m)p/V. Explicitly, 

\K{LI)JMj)= r'\afk\oi)\a{LI)JMj), (15) 
Jo 

with fK^{oi) a function given in terms of the Jacobi polynomials (see ap- 
pendix). The hyper-angular momentum K provides an upper bound for L 
and J,i.e.,J<L + I<K. Thus, if we take a value of K given by K^n = L^, 
we can argue that for K > tidal forces are less important than the forces 
between the fragments. On the contrary, when K < the forces between 
the fragments will be small compared to the tidal forces. The relative impor- 
tance of tidal forces compared to the forces between the fragments depends 
obviously on the values of Rat, Rbt and r. However, within the USA ap- 
proach, such relative importance between both types of forces is basically 
determined by the value of K. 

The Uncorrelated Scattering Approximation (USA) uses the expansion 
of the scattering wave function in terms of the hyper-angular momentum K. 
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Then, the Hamiltonian H is approximated by for the components of the 
wave function such that K > Km- Thus, tidal forces are ignored and the 
scattering is governed by the folding potential. This means that no excitation 
or break-up of the projectile occurs in these components. On the contrary, 
for K < Km the Hamiltonian H is approximated by H. Then, the correlation 
between the particles is ignored. It is very important to realize that the USA 
formulates different approximations for H in terms of the value of K, and 
not in terms of R. Thus, H (H^) is the approximate expression of H for any 
R, provided that K < Km {K > Km). 

It is important to formulate the USA to ensure that the interaction does 
not couple states with K < Km to states with K > Km- In order to do that, 
let P be the projector on the states with K < Km and Q the projector on 
the rest of states. The full Hamiltonians H can be expressed as HP + HQ. 
The USA implies that the term HP is approximated by PHP, while HQ is 
approximated by QH^ . This ensures that the time evolution of a state \i) 
is given by the sum of the evolution of P\i) and that oi Q\i), which remain 
mutually orthogonal. 

We can now study what our approximation implies regarding the S- 
matrix. We start with a regular solution of H and then we construct a wave 
packet. For t — > — oo the wave packet will be characterized by an incoming 
wave function times an internal state given by the ket \i) — \nLIJMj). The 
wave-packet for t +00 will be a product of outgoing waves times a combi- 
nation of states I/) = \n' L' I' J' M'j) , multiplied by certain coefficients, which 
are the matrix elements of the S-matrix operator between the states (/| and 

Then, 

{f\S\i)^{f\SP\i) + {f\SQ\i). (16) 

Within the USA model, the Hamiltonian H is replaced by H^ when referred 
to states with K > Km- This implies that the operator SQ can be approxi- 
mated by SpQ, where Sp is a c- number given by the elastic S-matrix for the 
calculation involving the folding potential. For K < Km the Hamiltonian 
H can be substituted for H, implying that the operator SP can be approx- 
imated by SP, where S is the S-matrix for the Hamiltonian H. As we will 
see in next section, S can be expressed in terms of the product of S-matrices 
of the two particles A and B, and it remains in the space of states with 
K < Km- Thus, SP ~ PSP. Finally, we can write 

{f\S\i)c^Sp5f,,+ {f\AS\i), (17) 
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where AjS" = P(S — Sf)P- Therefore, within the Uncorrelated Scattering 
Approximation, the S-matrix is given by the sum of two terms: the S-matrix 
coming from the folding model, which contributes only to the elastic scatter- 
ing, and a correction term that affects only to the components with K < Km, 
and which contains all the excitation and break-up effects. This term is given 
by the difference between the S-matrices of the two uncorrelated fragments 
and the S-matrix from the folding model. 



2.1 Boundary conditions 

Let us proceed now to describe the boundary conditions. Consider a regular 
solution of H^, characterized by a total energy E, orbital angular momentum 
L, internal angular momentum of the fragments /, total angular momentum 
J and internal energy e„: 

\KiljmAE)) = MPr,R)\n{LI)JMj). (18) 

The state \n{LI)JMj) can be written explicitly as 

\n{LI)JMj)^ p''dp(t>n{pMLI)JMj). (19) 



The channel wave function \n{LI)JMj) can be projected with the operator 
P, extracting the components with K < K^,- Thus, we have 

P\n{LI)JMj)= J2 {K\n)Li\K{LI)JMj), (20) 

K=L+I 

where the overlap is given by 

1/2 



{K\n)Li = / dpp 



da{p) 
dp 



fk'Hp)rMp)- (21) 



Thus, the asymptotic regular wave function, projected by P and written 
as an eigenstate of H^, becomes 

PKlijmAE)) = T.(K\n)u\KLijMjm (22) 

K 

where 

IKlijmAE)) = jK{Vn)\K{LI)JMj) . (23) 
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Using the Raynal-Revai transformation [2^, the incident wave function 
can be expressed in terms of the angular momenta associated to the coordi- 
nates -Rat; Rbt, 

\K{LI)JMj) = {LaLb\LI)jk\K{LaLb)JMj), (24) 

LaLb 

with J < La + Lb < K. Thus, we can write 

l^KLIJAlAE)) = Y: {LALB\LI).jK\^lL,LBJM,m ■ (25) 
LaLb 

Note that the Raynal-Revai coefficient vanishes for i^' < L + I or K < 
La + Lb- The state \^%claLbJMj{E)) ^ regular solution of for specific 
values of La, Lb and K 

\'^IlaLbJmAE)) = JK{Vn)\K{LALB)JMj) . (26) 

The angular momenta La and Lb are separately conserved in the scatter- 
ing process due to H. However, the hyper-angular momentum which 
is a good quantum number for is no longer conserved by H. Thus, 
one can proceed by using a new basis that keeps La and Lb as quantum 
numbers, but replaces K with other quantum number which is conserved 
by H. Note that this Hamiltonian ([12D conserves the energy, and hence 
the asymptotic momentum of each particle separately. In the appendix 
we show how to transform the states characterized by the values of K up 
to Km into states that have, approximately, a defined value of the mo- 
mentum of each particle. This transformation is achieved in terms of the 
Momentum Localized States (MLS) \1{LaLb)JMj) . These states depend 
on the momenta Pa and Pb which are strongly peaked around the values 
P| = V^MA/Mcos{/3e) and P| = V ^[m^Jm sm{/3e) , respectively. The en- 
ergies are given by Ej^ = {E — v) cos^(/?^) and E^ = {E — v) sin^i^Pi). The 
relation between the localized states and the original states is given by means 
of an orthogonal transformation, 

rii 

\K{LaLb)JMj) = Ym)LALB\^iLALB)JMj), (27) 
i=i 

where the number of momentum localized states, rii, coincides with the num- 
ber of states with definite K, rii = [{Km — La — Lb)/2] + 1. The coefficients 
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of the transformation are analytic expressions given in the appendix. Then, 
we can write 



l*M,JM,(^)) = Y.m)L,LML,L,JMAE)) . (28) 

1=1 

The state 1^°^^^ 

bJMj{E)) corresponds to a regular wave function in which 
the two particles A and B have linear momenta with narrow distributions 
around P4 and P^, and angular momenta La and Lb, respectively. If we 
define 

{£LALB\nLI)j= {K\n)u{LALB\LI)jKm)L^L,, (29) 

K=La+Lb 

which is a coefficient that depends on the bound wave functions and on 
analytic transformation coefficients, we can write finally, 

P\^^Slj,j^{E)) = Y: Y.{^LALB\nLI)j\^>',^^^^j^j^{E)). (30) 

LaLb «=1 

Note that this transformation relates the asymptotic states that define the 
boundary condition of the scattering problem, with the states for which each 
particle has a well defined angular and linear momentum. From this expres- 
sion, we can construct the incoming and outgoing waves just by making the 
adequate wave-packets. For t — > —00, eq. (^) relates the incoming parts 
of P\^fiLjMj{E)) and MlaLbJmAE)), while for t ^ +00, it relates the 
outgoing parts. 

2.2 S-matrix to bound and resonant break-up states 

Within the Uncorrelated Scattering Approximation, the two particles scat- 
ter independently inside the interaction region where the full Hamiltonian, 
projected on values K < Km, PHP, is replaced by PHP. The S-matrix is 
simply expressed in a basis of momentum localized states |^°L^LgjMj(-^))- 
A wave packet of these states at t ^ — 00 evolves according to PHP to give 
for t —>■ +00, the product of the S-matrices Sa{La, E^^)Sb{Lb, E^) times 
the wave-packet. Note that in writing this expression, one substitutes the 
narrow energy distributions of Ea and Eb of the MLS state for their central 
values 



10 



The matrix elements oi PSP in the channel basis \nLIJ) can be evaluated 
considering the transformation (|30|), 

{n'I'L'J\PSP\nILJ) = J2 Y.i^LALB\nLI)j{nL'I'\iLALB)j 

LaLb t=l 

X Sa{La,E'a)Sb{Lb,E's), (31) 

with J < La + Lb < Km- 

In order to evaluate the matrix elements for PSpP one should note that 
the operator Sp is a function of the orbital angular momentum L and the 
energy of relative motion En- It conserves the orbital angular momentum L 
and the internal angular momentum /, and is independent on K. Moreover, 
the operator P conserves L and / and projects on i^' < K^- Thus, we can 
write 

Km 

{n'l'L'J\PSFP\nILJ) = SriSviSpiL, E^) Y.{n'\K)iL{K\n)iL. (32) 

K 

that can be also expressed in the form, 

{n'I'L'J\PSFP\nILJ) = ^ Y^^^LaLb^LI) j{n' L' I'\ILaLb) jSf{L, E^), 

LaLb e=i 

(33) 

where we have used orthogonality properties. 

Then, the final expression for the S-matrix in the Uncorrelated Scattering 
Model results 

{n'I'L'J\S\nILJ) = dn'JndL'LSpiL, E) + {nI'L'J\AS\nILJ), (34) 

with 

rii 

{n'I'L'J\AS\nILJ)= ^ Y.'^^^aLb^LI) j{n' L' I'\ILaLb) j 

LaLb e=l 

X {Sa{La, E'^)Sb{Lb, E'b) - Sf{E, L)} . (35) 

This expression is valid for elastic scattering, inelastic scattering to bound 
states and break-up to resonant states in the continuum. We also notice that 
the expression (0) can be interpreted as the matrix element of the operator 
Sf + AS", with AS" = P{S — Sf)P, between the initial and final internal 
states. 
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2.3 Partial non-resonant breakup cross section 



In the previous section we have derived expressions for the S-matrix elements 
corresponding to final states, bound or resonant, that can be represented by 
normalizable wave functions. They are given in terms of the matrix elements 
of the operator AS* that is naturally described in the MLS basis \1{LaLb)J) , 
in which this operator is diagonal. Thus, only the overlap between the wave 
function of the final state and the states in the MLS basis is required to be 
known. The same procedure could be applied to calculate break-up to non- 
resonant continuum states. Provided that the corresponding wave functions 
are known in momentum representation, the overlap can be obtained. 

In this section we do not calculate the expressions for the break-up to spe- 
cific states in the continuum as they should depend on the detailed continuum 
states wave functions considered. Instead, we derive closed expressions for 
the non-resonant break-up cross sections, integrated over all the possible val- 
ues of the energies of the fragments, but characterized by a certain angular 
momentum of the fragments /' and of the relative motion L' . 

Making use of the completeness relation for the internal eigenstates, it is 
possible to derive a closed expression for the partial breakup cross section 
leading from the initial bound state \nLIJ) to all the final non-resonant 
continuum states characterized by the set of angular momenta {L', /', J}. We 
denote this cross section by a^{nLI —>■ LT). The details of the derivation 



are given in and will be published elsewhere. In this situation the final 



expression for the partial breakup cross section within the USA is given by 



afinLI^LT) = ^{2L + l){J2\{K{L'I')J\AS\nLIJ)f 

M) K 

-J2\{n'L'I'J\AS\nLIJ)\^}, (36) 

n' 

where Pq is the asymptotic incident momentum of the projectile. This expres- 
sion has a simple interpretation. The first term is the cross section induced by 
the operator AS to all the states labeled by the angular momenta K, L', I', J, 
that include the contribution of the bound and resonant states, which are 
explicitly subtracted by the second term. The summation with respect to K 
is extended to all the values between L' + /' and Km- 

It is also possible to obtain a compact expression for the breakup cross 
section corresponding to a total angular momentum J, a^finLI). This is 
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achieved upon summation of a^f{nLI LT) on the angular momenta 
L' and /' and taking into account the completeness property of the states 
\K[L' I')JMj) . This leads to the close expression 

(T^"(nL/) = ^{2L+l)[{nLIJ\{^SyAS\nLIJ)- ^ \{n' L' I' J\AS\nLIJ)\''}. 

M) n'L'V 

(37) 

Then, within the USA, the non-resonant breakup cross section for a given 
total angular momentum is calculated as the dispersion of the operator AS" 
in the ground state of the projectile, subtracting the contribution of the other 
bound and resonant states. 

3 Application to the d + ^^Ni reaction 

In this section we apply the uncorrelated scattering approximation to the 
analysis of elastic and breakup scattering of d by ^^Ni. Though the USA 
is expected to work better for more loosely bound projectiles such as or 
^^Be, for which the correlations between the fragments are weaker than for 
the deuteron, we start studying the case of the deuteron because this is a 
much better known system for which numerous calculations and experimental 
data already exist. 

The reaction d + ^^Ni has been extensively studied by the Kyushu group 
by means of Continuum Discretized Coupled Channel Calculations (CDCC) 
p9|- pl|] . It has also been used as a test case of the adiabatic approximation 
||32| and the Glauber multiple-scattering theory ||3^. All these approaches 
predict an important effect of the coupling to the breakup channels that 
results in a significant departure of their predictions compared to the folding 
model calculation. This is a characteristic phenomenon of reactions involving 
halo nuclei. Therefore, some of the conclusions arising from the analysis of 
reactions with deuterons can be also extended to the case of exotic nuclei. 

We first analyze the elastic scattering data at 80 MeV. As already men- 
tioned, the calculation of the elastic 5'-matrix elements within the USA re- 
quires the following ingredients: 

i) The internal wave function of the deuteron. Within the USA this wave 
function enters in both the folding potential and the coefficients {n\K)u ap- 
pearing in eq. (|34|) . We adopt in this work a simple model of the deuteron 
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which results from the assumption that the proton-neutron potential is sep- 



arable in momentum space [Q. In this model the S-wave component of 
the deuteron ground state is described in momentum space by the simple 
analytic expression 

where Uj^mC is related to the range of the proton-neutron interaction, B is 
the binding energy of the deuteron [B = 2.22 MeV) and is a normalization 
constant. We neglect the small D-wave component of the ground state wave 
function and the proton and neutron intrinsic spins. 

ii) The second ingredient of the USA refers to the two-body S matrices for 
the constituents. In the case of the deuteron, the proton-target and neutron- 
target S'-matrices {Sp, Sn) are required for values of the angular momenta 
in the interval < Lp, < K^, and values of the energies determined by 
the momentum localized states, which lie in the range < Ep, < E — v. 
These S matrices have been calculated by means of optical potentials as it is 
done in the CDCC, adiabatic and Glauber calculations. In these formalisms 
the proton and neutron optical potentials are evaluated at half of the inci- 
dent deuteron energy and so, the energy dependence of the optical potential 
parameters is neglected. This approximation is based on the assumption 
that the proton and neutron move approximately with the same velocity of 
the deuteron center of mass, and the dispersion around this value is small. 
Those configurations for which one of the fragments carry the whole available 
energy must be highly suppressed. Within the USA this fact is explicitly in- 
cluded in the coefficients, (n(L/) J|£L^Lb), that can be physically regarded 
as the amplitude probability of having the constituents of the projectile with 
angular momenta La and Lb, and energies E^ and E^ within a state char- 
acterized by the relative angular momentum L, internal momentum / and 
total angular momentum J. In fact, these coefficients favour those configu- 
rations for which each one of the particles carries half of the incident angular 
momentum and half of the available energy, according to a classical picture. 
It is important to note that, within the USA, the energy dependence of the 
optical potential is naturally taken into account by solving the proton-target 
and neutron-target Schrodinger equations at definite scattering energies. In 
particular, we perform the calculations with the optical potentials of Ref. ||35|| . 

The USA also requires the introduction of two parameters, the average 
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potential, v, and the cut-off hyper angular momentum, Km- These param- 
eters have a clear physical interpretation and so they could be set to some 
reasonable values, or otherwise fitted to improve the agreement with exact 
calculations. Nevertheless, we keep in mind that our purpose in developing 
the model is to apply it to weakly bound systems for which K.^ must be large 
and V small. Then, we have adopted the simplest approach. We take 
large enough to achieve convergence for the S-matrix elements. This means 
that the correlations between the fragments can be neglected even at large 
distances. For the average potential we take w=0. Several test calculations 
have revealed a weak dependence on this parameter. Hence its choice does 
not affect significantly the results. 

The calculated angular distribution of the elastic differential cross sec- 
tion, divided by the Rutherford cross section, is plotted in Fig. |l|. Experi- 
mental data are represented by circles. The dotted line refers to the folding 
calculation, which clearly overestimates the cross section at intermediate an- 



gles. The dashed line is the result of the CDCC calculation taken from ||29 
that, compared with the folding model, predicts a significant reduction of the 
cross section at intermediate angles. This is a consequence of the coupling 
to breakup channels. The solid line corresponds to the USA calculation with 
the cutoff hyper angular momentum ii'm=30, for which convergence of the S- 
matrix elements is achieved. As shown, experimental data are not accurately 
reproduced by the full USA calculation that overestimates the reduction of 
the cross section with respect to the folding model due to break-up effects. 

In order to provide an explanation of this result we investigate the partial 
breakup cross section, Uj". Within the USA, this quantity can be easily eval- 
uated without any explicit description of the continuum states, by means of 
eq. (|37D that only requires the introduction of the ground state wave function. 
In Fig. |2| we show the distribution of cTj" versus the total angular momen- 
tum J calculated in the USA (thin solid line). This result is compared with 
the CDCC calculation (dashed line) performed within the subspace /'=0, 2. 
Thus, in order to enable a meaningful comparison between both approaches, 
we have also calculated within USA, the contribution to the breakup cross 
section due to the S (/'=0) and D (/'=2) components. The result, plot- 
ted in Fig. PI by the thick solid line, shows a fairly good agreement with 



the CDCC calculation (dashed line) p9|,^. Apart from the agreement in 



the overall magnitude, the angular momentum dependence is also accurately 
reproduced, including the surface-peak nature of the elastic breakup process. 
From the results for the partial breakup cross section (Fig. 2), it is clear 
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Figure 1: Elastic differential cross sections angular distributions (as ratio 
to Rutherford) for d + ^^Ni scattering at 80 MeV. The dotted, dashed and 
thin solid line correspond, respectively, to the folding, CDCC and full USA 
calculations. The thick solid line is the restricted USA calculation which 
contains only the effect of /' = 0, 2 breakup states. Experimental data 
are given by circles. 
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that the USA predicts an important contribution to the breakup coming from 



continuum states with internal angular momenta I' > 2. Austern et al. [ST 
have studied the convergence of the CDCC calculation for this reaction with 
respect to the cut-off internal angular momentum. Their calculations reveal a 
small contribution to the breakup cross section coming from breakup channels 
with /' > 2. For example, within the model space /'=0,1,2,4,6, for which 
a good convergence of the solution is achieved, the breakup cross section 
associated to J=17 is 17.27 mb, whereas the USA predicts a value of 22 mb. 
We interpret this discrepancy as a consequence of the basic approximation 
involved in the USA, i.e., to neglect the correlations between the constituents 
in the scattering process. In this sense, our treatment is opposite to the 
adiabatic approximation, in which the internal coordinate is assumed to be 
frozen during the collision, keeping the constituents strongly correlated, and 
thus avoiding the breakup to high angular momentum states. 

This effect could also explain the discrepancy encountered for the differ- 
ential elastic cross section. In order to provide a numerical assessment of this 
hypothesis we have performed a new calculation for the elastic cross section 
in which the contribution of the continuum channels with J' 7^ 0, 2 has been 
excluded in an effective way. We recall that the elastic S-matrix in the USA 
is given by the sum of two terms, the first one coming from the folding po- 
tential and the latter, AS", which describes dynamic polarization effects due 
to the coupling to breakup channels. Therefore, this second term arises from 
the effect of the tidal forces. However, the tidal potential vt has vanishing 
diagonal matrix elements on the ground state of the projectile. Thus, the 
contribution of break-up states with angular momenta /' to the elastic matrix 
elements of AS' (up to lowest order in the tidal forces) depends on second 
order coupling through the expression 

{(j)oLIJ\AS{I')\(f)oLIJ) oc {(poLIJ\vT\bu;L'I'J)G+{bu;L'I'J){bu;L'I'J\vT\(j)oLIJ), 

bu,L' 

(39) 

where G'^ {bu; L' I' J) is a propagator and the sum extends to the breakup 
states. In the particular case in which off-shell dependence is removed from 
the the breakup states involved in the propagator, the right hand side of 
eq. (^9|) is proportional to the square of the distorted wave integral of the 
tidal potential. To lowest order, this term is simply proportional to the 
break-up cross section. 

A restriction in the number of breakup states considered will produce 
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a reduction in the elastic matrix elements of AS. The results presented 
in Fig. 2 show that the USA calculation including /' = 0,2 break-up is 
correct, whereas the calculation of break-up to larger angular momentum 
is overestimated. Thus, making use of the proportionality between elastic 
matrix elements of AS' and the break-up cross section, we get 

bu / T' — 2) 

{nLIJ\AS{r = 0, 2)\nLIJ) = {nLIJ\AS\nLIJ) ^ ^ ' ^ (40) 

where crj"(/' = 0, 2) is the restricted breakup cross section and AS{I' = 0, 2) 
its associated elastic S-matrix contribution. The corresponding elastic differ- 
ential cross section resulting from this prescription, which we call "restricted" 
USA calculation, is represented by the thick solid line in Fig. |1|. The result 
is in very good agreement with the experimental data, supporting the hy- 
pothesis that the strong coupling to high spin breakup states within USA 
is responsible for the reduction in the elastic cross sections. If this coupling 
to high spin states (/' > 2) is excluded from the calculation, and the elastic 
matrix elements are modified accordingly, then both the elastic differential 
cross sections and the partial break-up are well described. 

An interesting question is to assess the validity of the USA at low energies 
where other models, which are successfully applied to the high energy regime, 
fail to reproduce the experimental data. This is the case of the sudden 
approximation and, in particular the Glauber model, that can be regarded 
as a high energy approximation to the adiabatic treatment. To this end we 
have applied the USA to the reaction d + ^^Ni at 21.6 and 56 MeV, for which 



experimental data are available ||33|. We use the same wave function for 



the deuteron ground state as in the previous calculations. The Percy optical 
model parameterization |^ for protons and neutrons has been selected. The 
angular distribution of the elastic differential cross section is shown in Figs. 
^ and ^ where we compare the experimental data with the USA result (thin 
solid line), the USA result considering break-up with J' = 0, 2 (thick solid 
line), the CDCC calculation (dashed line) and the folding model (dotted 
line). 

The results obtained for 56 MeV, shown in Fig. are qualitatively similar 
to the ones for 80 MeV (Fig. 1). The full USA calculation displays a too 
large reduction in the elastic cross section with respect to the folding model 
prediction. However, the "restricted" USA calculation, which considers only 
the effect of break-up to I' = 0,2, reproduces accurately the experimental 
data. 
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J 

Figure 2: Partial breakup cross section, as a function of the angular mo- 
mentum. The thin solid line is total breakup cross section predicted by the 
USA calculation (eq. The thick solid line corresponds also to the USA 
calculation, but including only the S and D continuum channels. The dashed 
line is the analogous calculation in the CDCC approach. 
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Figure 3: Elastic differential cross sections angular distributions for d + ^^Ni 
at 56 MeV. The meaning of the curves is the same as in Fig. |l|. 



In Fig. H we present the results for 21.6 MeV. Here there is a fairly good 
agreement between the full USA calculation and the experiment up to 100 
degrees. It is noticeable that the angular region between 50 and 100 degrees 
is even better described by the USA model than by the CDCC approach. The 
"restricted" USA calculation improves the agreement with the experimental 
data. 

These results indicate that the USA is adequate to calculate deuteron 
break-up cross sections to states with /' = 0, 2, which are the most important 
break-up components, although it overestimates the break-up cross sections 
to deuteron states with /' > 2. As far as elastic scattering is concerned, the 
full USA calculation gives in general, a too strong reduction in the cross sec- 
tions. However, the "restricted" USA calculation, which takes effectively into 
account only the effect of coupling to /' = 0, 2 break-up states, reproduces 
satisfactorily deuteron elastic scattering data at all the energies considered. 

Although more detailed tests should be required in order to delimit the 
range of validity of the USA, these preliminary results suggest that the model 
can be used as an alternative tool to analyze experimental data of halo nuclei 
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at relatively low energies, around and above the Coulomb barrier. 

4 Summary and conclusions 

We have developed a new approach, called the Uncorrelated Scattering Ap- 
proximation (USA), applicable to the scattering of weakly bound nuclei. Our 
first requirement is that the target is very heavy compared with the fragments 
of the projectile. Although this assumption is not essential to the model it 
simplifies importantly the general formalism. The USA arises from the fact 
that the three-body Hamiltonian corresponding to a composite projectile can 
be written as a sum of two non-interacting two-body Hamiltonians in two 
limit cases: i) when the tidal forces are neglected, and ii) when the forces 
between the fragments of the projectile are neglected. In the first case, the 
projectile remains in its ground state, the interaction with the target is given 
by the folding potential and the S-matrix is given by the solution of the cor- 
responding two-body scattering problem Sf{L,E). In the second case, the 
particles scatter independently. This means that if the initial wave function 
is given by the product of incident waves of the particles, the wave function 
in the final state is characterized by a product of outgoing waves, multiplied 
by the corresponding S-matrices, Sa{La, Ea)Sb{Lbj Eb). Thus, within the 
USA one neglects tidal forces in the scattering region corresponding to large 
separations, whereas in the scattering region, corresponding to small sepa- 
rations, the forces between the fragments of the projectile are ignored. The 
regions of large and small separations are defined in terms of the hyper an- 
gular momentum K. Large separations correspond to i^' > K„i, and small 
separations to K < Km, with Km defined such that its turning point corre- 
sponds to a distance Rm for which tidal forces and the forces between the 
fragments are comparable. For K > Km the forces between the fragments 
dominate, while for K < Km tidal forces are more important. Thus, Km is 
a parameter of the USA model. The other parameter, v, is a constant that 
substitutes the interaction between the fragments when K < Km- 

The proper boundary conditions for the region K < K^ requires the 
projection of the initial state in a discrete basis of states which are strongly 
localized for certain values of the relative momentum of the fragments. These 
momentum localized states can be transformed analytically to states which 
have definite values of the hyper-angular momentum K. Then, the Raynal- 
Revai coefiicients are used to transform the initial states, given in terms of the 
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projectile-target angular momentum L and the internal angular momentum 
/, to states that depend directly on the angular momentum of the fragments 
with the target, La and Lb- Finally, a certain combination of i^'- values is 
done in order to obtain states with strongly localized values of the kinetic 
energy of each fragment. This procedure allows to write the S-matrix of the 
projectile in terms of the product Sa{La, Ea)Sb{Lb, Eb) of the S-matrices of 
the fragments, evaluated at certain values of the energy and angular momen- 
tum of the fragments with respect to the target. Then, the formalism can 
be applied directly to describe scattering to bound and resonant break-up 
states, as well as to direct break-up states. 

The main advantage of the USA is that it allows to express the S-matrix 
of the composite system in terms of the wave function of the bound states 
and the scattering S-matrices of the fragments with the target. These S- 
matrices can be easily obtained from calculations that fit the cross sections 
of scattering experiments of the fragments with the target. In this sense, the 
fragment-target interaction is treated to all orders, and the only assumption 
one makes on the reaction dynamics is to neglect the correlation between the 
fragments. Thus, the USA approach presents some advantages compared 
with previous treatments. Effects of absorption or "shadowing", which are 
included by means of ad— hoc profile functions in some spectator models |||,||, 
appear naturally here, as occurs in Glauber approaches ||22|,^, because the 
S-matrices Sa{La,Ea) and Sb{Lb, Eb) of the fragments include the effects 
of absorption for low La and Lb values. Compared with semi-classical ap- 
proaches p4|,|25|, the USA does not make use of the concept of classical tra- 



jectories for the relative motion. Compared with Glauber analyses ||2^,|23 



the USA justifies the expression of the three-body S-matrix in terms of the 
product of the S-matrices of the fragments from very general considerations, 
that do not require to make use of the eikonal assumption of forward scat- 
tering, or straight line trajectories. Also, our expressions depend on the 
S-matrices of the fragments evaluated at integer orbital angular momenta, 
while previous Glauber-type approaches rely on evaluating by interpolation 
the S-matrices of the fragments as a function of real impact parameters, which 
would correspond to non-integer L-values. Finally, the USA also differs from 
Glauber-type approaches, and other approaches based on the adiabatic ap- 
proximation, for which the velocity of all the fragments with respect to the 
target is equal to the projectile velocity, and so the energies of the fragments 
are fixed quantities proportional to the masses. In the USA, the relative 
energy of the fragments with respect to the target, Ea, Eq can take different 
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values, reflecting the fact that in the incident projectile the fragments have 
a certain momentum distribution. As the USA does not make use of the 
eikonal approximation, which is essentially a forward-scattering approach, 
we consider that it can be useful when applied to scattering of weakly bound 
halo nuclei at low scattering energies, comparable to the Coulomb barrier. 

Apart from the theoretical formulation of the USA, we have also presented 
a preliminary application of the approach to the reaction d \ ^^y.i at 21.6 
and 80 MeV. The parameter Km is set to a large value, and v is taken as 
zero. This corresponds to neglect completely the effect of proton-neutron 
correlation in the scattering. Concerning elastic scattering at 80 MeV, the 
USA gives rise to a depletion of the cross sections at intermediate angles 
with respect to the folding model prediction. The experimental data and 
the CDCC calculations also show this depletion, although not so large as 
in the USA model. The break-up cross sections at 80 MeV to /' = and 
I' — 2 states in the USA approach is consistent with the CDCC calculation. 
However, break-up cross sections to /' > 2 states, which are sizeable in the 
USA, are strongly suppressed in the CDCC calculations. When the USA 
calculation of the elastic scattering is modified to exclude effectively the 
coupling to /' > 2 states, the experimental differential cross sections and 
CDCC calculations are very well reproduced. Using the same procedure for 
the elastic cross sections at 21.6 and 56 MeV, the experimental data are also 
accurately reproduced. 

Our interpretation of these results for deuteron scattering is that the 
proton-neutron correlations, which are neglected in the USA calculation, play 
a significant role. These correlations show up in the evaluation of break-up 
cross sections to /' > 2 states. Within the USA model, the proton and 
neutron scatter independently from the target, and so, after the scattering 
they have a certain probability to end up in a state with large relative angular 
momentum. However, in reality the proton and neutron remain correlated, 
and this correlation suppresses large /' break-up components in the wave 
function. Once these components are excluded, both elastic and break-up 
cross sections are well reproduced by the USA calculation. 

We expect that, for the scattering of halo nuclei at energies around the 
Coulomb barrier, the effect of the correlations between the fragments will be 
less important than in the case of deuteron. Thus, the uncorrelated USA 
calculation will be more reliable and we may expect to see more clearly the 
features of the USA calculations. These include an important decrease of 
the elastic cross section with respect to the folding model prediction, and 
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sizeable contributions to break-up with large angular momentum between 
the fragments. 
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Appendix A 

Momentum Localized States associated to hyper-spherical 
harmonics (HH) 

In this appendix we show how states characterized by definite values of the 
linear momenta of the constituents can be built starting from a truncated 
basis of hyper-spherical harmonics. 

First, let us introduce the HH states, \K{LaLb)JMj) , characterized by 
the hyper-angular momentum K, the orbital angular momenta La and Lb 
of particles A and B, respectively, and the total angular momentum J and 
its projection Mj. We define the angle f3 that connects the linear momenta 
of the particles with the hyper-momentum V through the relations Pa — 



I MaIM)V cos (5 and Pb = (y/ Mb /M)V sin (3. The states \K(LaLb)JMj) 
can be expressed in terms of states with definite values of /? as 

\K{LaLb)JMj) = r^\(3fji^''''{f3)\(3{LALB)JMj), (A.l) 
Jo 



where 



= A^-^^^(cos/5)^-^+^(sin/3)^^+^pi^^+^'^^+^\cos2/3), (A.2) 

where P^"'^^ denotes a Jacobi polynomial of degree n — {K — La — Lb)/2 
and N^^^'^ represents the normalization constant 



LaLb 



2n\{K + 2){n + LA + LB + l)\ 
V{n + La + |)r(n + Lb + 1) 



(A.3) 
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The function is normalized to unity in the variable in the vari- 

able /?: 

^^\(3[fk^'^-{(3)f = 1. (A.4) 
L J 

This property also guarantees that the HH are normalized to unity. The 
hyper-angular momentum can take the values K = La + Lb, La + Lb + 
2,... giving rise to an infinite number of states in eq. (|A.1|) . One of the 
basic ingredients in the Uncorrelated Scattering Approximation developed in 
Section III is the introduction of a maximum hyper-angular momentum Km 
which reduces the infinite set of states to a finite one: K = La + Lb, La + 
LB + 2,...,Km- In terms of the index n we have the subset: n = 0, 1, A^ — 1, 
with A^ = {Km -La- Lb)/2 + 1. 

The family of Jacobi polynomials, {P^"'^^(x); n = 0, ...,A^ — 1}, constitute 
an orthogonal set of functions in the interval (—1, +1) with respect to the 
weight function uj{x) = H^Ti^Y: 

^\xu{x)[Pi'''''\x)]^ = hr, (A.5) 
1 

where 

^ T[n + a + l]T[n + b+l] 

" 2{2n + a + b + l)n\T[n + a + b+iy ^ ' ' 

From this family of Jacobi polynomials a new set of A^ polynomials of 
degree A^ — 1 are defined as p9|] 



Km 

E 

K=La+Lb 



AT, 



LaLb 



K 



P^^)(x)P^^)(x,); i = l,...,N. (A.7) 



Here, {x^; (i = 1,2,..., A^} represent the zeros of the polynomial Pj 



N ' 



X] 



Direct application of Christoffel-Darboux formula (see ref. ||4^) on the pre- 
vious expression leads to 



N 



X) 



{a,b) . 



N 



(X — , 



(A.8) 



where kN is the coefficient of x^ in P^''''{x). The polynomials Q^l'!^^(x£, x) 
are orthogonal in the same interval and relative to the same weight function 
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as the original polynomials. The expression above shows that the polynomial 
QN-ii^e^ ^) vanishes at the points x = Xg] s = 1, . . . , N except at x = X£. 

The new set of polynomials allows us to construct Momentum Localized 
States (MLS) which are defined as 



\i{LALB)JMj) = r^'rf/3(cos/3)^^+i(sin/3)^^+i 

Jo 

X Q^^^^'''''^'''^\x,,cos{2P))\P{LaLb)JMj) (A.9) 



with 

n -1 



(LaLb) 



(A.IO) 



Substituting the explicit expressions of /itv-i, ^at-i and [0 (|A.10|) 
one gets 

(LaLs) _ 1 {2N + a + bY T[N + a + l]T[N + b+l] 1 - xj 



23 (AT + a)2(iV + 6)2 mT[N + a + b+l] 



2 ' 



(A.ll) 

where a = Lb + 1/2 and b = La + 1/2 in our case. 

The state \1{LaLb)JMj) is characterized by values of the momenta of 
the two particles sharply peaked around P'x = {\J Ma/ M)V \J (1 + x^)/2 and 

Pb = {\J Mb/M)V\J (1 — X()/2. This localization is enhanced as the num- 
ber of states increases which, in turn, depends on the cut off hyper-angular 
momentum K^. 

An orthogonal transformation can be defined between the truncated hyper- 
spherical basis and the MLS basis: 

\e{LALB)JMj) = {K\e)LALB\K{LALB)JMj). (A.12) 

K=La+Lb 

This relation can be inverted allowing to express the HH in terms of the 
MLS: 

N 

\K{LaLb)JMj) = J2{i\K)LALB\^iLALB)JMj). (A.13) 

1=1 
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The coefficients of the transformation can be computed from the definition 
of the MLS and the polynomials QN-i{xe,cos{2l3)): 



imL^LB = ")'/'iV^-^''"Pf (A.14) 

In Fig. H the functions /^'*^^(/?) (upper part) and the corresponding 
MLS (lower part) are plotted versus the variable /? for the case La=Lb=5. 
A basis with A^=5 states has been used, corresponding to the values of the 
hyper-angular momentum 7^^=10 to 18, in units of two. The MLS are labeled 
with the index i which is associated with the roots of the Jacobi polynomial 

(5-1-1^5+1) 

P5 ^' ^ (x). Thus, the MLS £ = 1 is localized around x = —0.64, that cor- 
responds to P = 1.13. In terms of the energy of the two particles this means 
that the fraction of the total available kinetic energy carried by particles A 
and B are ^a/(^ - v) = 0.18 and Eb/{E - v) = 0.82, respectively. This is 
in fact the most asymmetric situation for this value of Km- On other side, 
the MLS i = 3 corresponds to a physical situation in which both particles 
carries half of the total kinetic energy. 

In a similar way, it is possible to define HH with definite values of the 
hyper-angular momentum K, orbital momentum L, intrinsic spin I and to- 
tal angular momentum J and projection Mj. These new HH, denoted by 
\K{LI)JMj), can be related to the set \K{LaLa)JMj) by means of the 



Raynal-Revai transformation ||27 



\K{LI)JMj) = {LaLb\LI)jk\K{LaLb)JMj). (A.15) 

LaLb 

Introducing the angle a, which denotes the ratio between the internal 



momentum p and the hyper-momentum V, p = Pym/M sin a, it is possible 
to expand the new HH as 

\K{LI)JMj)= r^^dafk'{a)\a{LI)JMj), (A.16) 
Jo 

where the function /^^(a) is given by the analogous to (|A.2| ), i.e. 



f^\a) = iV^'(cosa)^+i(sina)^+ipy^^'''+^\cos2a), (A.17) 

with n = {K — L — I)/2 in this case. 

Introducing a cut-ff hyper-angular momentum, K^., it is possible to con- 
struct a set of Momentum Localized States associated with the truncated 
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Figure 5: /^^^''(/S) versus the variable /3 for La = Lb=5 and ir=10,12, 14,16 
and 18 (upper figure) and associated Momentum Localized States, labeled 
by the index i, defined in the text (bottom figure). 
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basis \K{LI)JMj) , L + I < K < K^^. They are developed in terms of the 
angle a as 

\j{LI)JMj) = {w^'^Y^/^^ da{cosa)''+\smay+^ 

X QN-!'^^^\xj,cos2a)\a{LI)JMj) (A.18) 

where N is the number of values of K of the HH truncated basis and {xj;j = 

1,...,A^} represent the zeros of the polynomial Pj^ ^' ^ (x). The state 
\j{LI)JMj) is characterized by narrow distributions of the relative and center 
of mass momenta around the central values pj — \f^'P^{l — Xj)/2 and 

= V^{l + x^)/2. 

The original basis and the MLS basis are related by means of the orthog- 
onal transformation 

N 

\K{LI)JMj) = J2{j\K)li\j{LI)JMj), (A.19) 

where the transformation coefficients are given by 

imu = {w'N'j)'^'Nk'Pi'^'^'^'^'^'\xj). (A.20) 
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